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ABSTRACT 

Hierarchical structure formation theory is based on the notion that mergers drive galaxy evolution, 
so a considerable framework of semi-analytic models and N-body simulations has been constructed 
to calculate how mergers transform a growing galaxy. However, galaxy mergers are only one type of 
major dynamical interaction between halos - another class of encounter, a close flyby, has been largely 
ignored. We use cosmological N-body simulations to reconstruct the entire dynamical interaction 
history of dark matter halos. We present a careful method of identifying and tracking a dark matter 
halo which resolves the typical classes of anomalies that occur in N-body data. This technique allows 
us to robustly follow halos and several hierarchical levels of subhalos as they grow, dissolve, merge, 
and flyby one another - thereby constructing both a census of the dynamical interactions in a volume 
and an archive of the dynamical evolution of an individual halo. In addition to a census of mergers, 
our tool characterizes the frequency of close flyby interactions in the Universe. We find that the 
number of close flyby interactions is comparable to, or even surpasses, the number of mergers for 
halo masses > 10^^ Mo/h at z < 2. Halo flybys occur so frequently to high mass halos that they are 
continually perturbed, unable to reach a dynamical equilibrium. In particular, we find that Milky 
Way type halos undergo a similar number of flybys as mergers irrespective of mass-ratio for 2; < 2. 
We also find tentative evidence that at high redshift, z > 14, flybys are as frequent as mergers. Our 
results suggest that close halo flybys can play an important role in the evolution of the earliest dark 
matter halos and their galaxies, and can still influence galaxy evolution at the present epoch. Our 
companion paper quantifies the effect of close flyby interactions on galaxies and their dark matter 
hosts. 

Subject headings: cosmology: theory — cosmology: dark matter — cosmology: large-scale structure 
of universe — galaxies: evolution — galaxies: halos — galaxies: interactions — 
methods: numerical 



1. INTRODUCTION 

In a ACDM Universe, the smallest dark matter 
halos form first; bigger halos are then formed via 
successive mergers with smaller halos. Thus, mergers 
are instrumental in the formation and evolution of 
halos. Galaxy mergers are rare, punctuated events 
in a galaxy lifespan that nonetheless dr amatically 
change it - from its mo rpholo gy fe.g., iHolmberd 
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2007; Lin ct al. 2008; 'T acconi et al] 120081: lUrrutia et a!,. 
200 8; R yan ct a l. .2008 [k Collisionless cosmological 
N-body simulations can be used to measure halo merger 
rates, where a merger is defined to occur when a bound 
dark matter halo falls into another bound dark matter 
halo. Various simulations measuring such halo merger 
rates agree to within a factor of 2 (see Hopkins et al] 
I20101 for a discussion on merger rates from collisionless 
simulations). Galaxy merger rates can then be inferred 
from the subhalo mergers within a pri mary halo by 
assum ing a Mha,in — M^^ ] relation (e.g., [Guo fc WhTta 
[200l I Wetzel et all [200a iBehroozi et all 120101 ) or di- 
rectly measured in hydrodyn amic simulations (e.g., 
iMaller et all 120061: iSimha et all [2009.) . Observationally, 
merger rates are typically derived from close-pair counts 
- i.e. galaxies with small projected separations and 
relative velocities - and are globalized using an estimate 
of the lifetime or du ration of the observed merger 
phase (iLotz et al.ll2010[ ). 

Ultimately, galaxy mergers are successful in shaping 
galaxy properties because they cause a large perturba- 
tion within the potential. Even an orbiting satellite can 
distort the underlying smooth galaxy potential and pro- 
duce observable effects. For example, the HI warp in the 
Milky Way disk m ay be tidally triggered by the Large 
Magellanic Cloud (iWeinberg fc Blitz ,2006) . However, 
one entire class of galaxy interactions also capable of 
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causing such perturbations - galaxy flybys - has been 
largely ignored. 

Unlike galaxy mergers where two galaxies combine into 
one remnant, flybys occur when two independent galaxy 
halos interpenetrate but detach at a later time; this can 
generate a rapid and large perturbation in each galaxy. 
We developed and tested a method to identify mergers 
and flybys between dark matter halos in cosmological 
simulations and to construct a full 'interaction network' 
that assesses the past interaction history of any given 
halo. With this new tool, we are undertaking the first 
systematic study to quantify the frequency of flybys and 
its effect on galaxy evolution. 

In this paper, we present our technique for determin- 
ing the network of dynamical interactions for halos in an 
N-body simulation, and present a census of halo flybys 
and mergers in the Universe. We discuss our simulation 
and halo- finding technique in Section 2. Section 3 de- 
scribes our technique to identify halo flyby interactions, 
and to construct a halo interaction network. Section 4 
presents the results, and Section 5 discusses the implica- 
tions of flyby interactions and previews the next paper 
in the scries. The appendix covers our method of link- 
ing parent and child halos, including ways to mitigate 
common problems that plague this process. 

2. SIMULATION TESTBED AND HALO IDENTIFICATION 

We use a high-resolution, dark matter simulation with 
1024^ particles i n a box of length 50 Mpc/h with WMAP- 
5 cosmological (jKomatsu et al.l l2009) parameters as a 
testbed to develop our technique. The initial particle 
distribution is obtained from a Zeldovich linear approx- 
imation at a starting redshift of z = 249 and is evolved 
using the adaptive tree-code, GADGET-2 (iSprin gel et al.l 
HjDOlb; Springcl 2005). The dark matter particles have 
a fixed gravitational co- moving softening length of 2.5 
kpc/h. We store 105 snapshots spaced logarithmically 
in scale-factor, a — 1/(1 -I- z), from z = 20 to 0. This 
translates into a timing resolution < 50 Myr for z ^ 3 
and ~ 150 Myr for z < 3. Since the fundamental mode 
goes non-linear at z = 0, we will only present results up 
to z = 1 where the 50 Mpc/h box is still a representative 
cosmological volume. 

In addition, we ran two simulations designed to explore 
the effect of mass and timing resolution on the fiyby phe- 
nomenon. These two volumes, with 512^ and 1024^^ par- 
ticles respectively, have the same size and cosmology as 
our testbed simulation and were evolved from z — 249—1. 
The particles were drawn from identical phases, and with 
161 snapshots total, had timing resolution better than 
200 Myrs. 

To begin identifying halos, we first use a Friends-of- 
Friends(FOF) technique with a cano nical linking length 
b = 0.2 (~ 10 kpc/h) pavis et al.l [l985). We require 
at least 20 particles (~ lO^MQ/h) to define a halo, but 
our halo interaction network uses only those halos with 
greater than 10 pa rticles; we discuss the effects of this 
limit in Section [A. 21 Subhalos (down to multiple hierar- 
chy levels) are identifi ed using the SUBFIND algorithm 
(jSpringel et al.ll2001a[ ). The SUBFIND algorithm identi- 
fies subhalos as bound structures around a density max- 
im 43. The remaining particles in the FOF halo are com- 
prised of particles that are bound to the potential of the 
main halo but not bound to any subhalo. For the remain- 



der of this paper, all references to the main or primary 
halo will mean this bound set of "background" parti- 
cles. For deta ils on the SUBFIND algorithm, we refer the 
reader to Spr ingel et al.l (|2001al ). We caution that SUB- 
FIND, or any density-based halo extraction technique, 
only recovers subhalo masses enclosed within a region of 
higher subhalo densi ty compared to the background (see 
iMuldrew et al.|[20To( ) . Therefore, a subhalo will lose mass 
depending on the radial distance from the center of the 
main halo; i.e., the subhalo mass is only to be trusted 
modulo a density contrast. If we want a detailed look 
at the subhalo mass loss evolution, then we must care- 
fully correct for this bias. Fortunately, at this stage we 
simply require the existence of the subhalo; we caution 
that there is one consequence of this bias that presents 
a challenge: in the most extreme case, a subhalo pass- 
ing through the center of the main halo can disappear 
for multiple snapshots only to reappear later, where it 
can be mistaken for a new subhalo. We account for such 
missing subhalos while constructing the halo interaction 
network (see Section ETTj) . Since our focus is on flybys, 
we relegate most of the discussion of our technique to the 
Appendix, including resolution tests and comparison to 
other techniques. 

3. TAXONOMY OF HALO DYNAMICS 

In a collisionless simulation, dark matter halos can 
grow via mergers or through smooth accretion. They 
can also lose mass through dynamical interactions, such 
as flybys or 3-bod y encounters (e.g. iSales et all 120071 : 
iLudlow et all[200l : primary halos can even be stripped 
of entire subhalos. Halos close to the numerical reso- 
lution limit can disappear for multiple snapshots. A 
primary halo can fall into another primary halo and 
continue to survive as a subhalo that is subsequently 
stripped. Ideally, a robust account of a halo's dynamical 
past captures all such scenarios and produces a physically 
meaningful history for individual halos. In a nutshell, our 
halo interaction network attempts to capture all of these 
interactions. We construct our halo interaction network 
in a three-step fashion: 1) for every halo, we flnd the 
best child halo at some later snapshot (if it exists); 2) we 
trace the complete set of parent-children pairs through 
all snapshots to construct the halo network; and 3) us- 
ing the phase space history, we characterize the type of 
interaction for each pair. 

Capturing the rich history of interactions of dark mat- 
ter halos requires identifying every possible interaction. 
Most of the interesting phenomena begin with a main 
halo falling into another main halo. Assuming that a 
secondary main halo falls into a primary main halo, the 
following outcomes are possible: 

• The secondary (now a subhalo) continues to orbit 
inside the main halo and ultimately disrupts inside 
that same main halo (merger started). 

• The secondary continues to orbit the primary halo 
and exists at z = (merger completed). 

• The secondary merges with another subhalo inside 
the same primary halo (subhalo merger). 

A particle is associated with only one subhalo 
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• The primary main becomes a subhalo itself by 
merging with another halo. The secondary is now 
a sub-subhalo (merger). 

• The secondary halo enters and then detaches from 
the primary halo (grazing flyby). 

• Either the secondary or the primary disappear 
(transients) 

It is also possible that a subhalo appears without ever 
having been a main halo. Based on our snapshot timing 
resolution work (see Appendix), we find that these sub- 
halos formed as main halos very close to the primary and 
merge in between snapshots. 

An interaction begins the moment a halo changes its 
state, e.g., when a main halo becomes a subhalo of an- 
other main halo. We categorize each interacting pair of 
halos by tracking the behavior of the pair into the future. 
At the moment the encounter begins, we tag the interac- 
tion according to the encounter type. For instance, when 
a subhalo enters a primary and eventually disrupts inside 
it, the interaction is tagged as a merger. Along with this 
tag, we store the redshift of infall (when the subhalo first 
becomes part of the main) and the destruction redshift. 
Thus, we automatically get a duration for every inter- 
action. In addition, we store the minimum separation 
between the subhalo and primary halo and the normal- 
ized radius at destruction - this allows us to take census 
of the galactocentric distances reached by any particular 
interaction as a function of primary mass, mass ratio and 
infall redshift while accounting for the dynamic behavior 
of the primary halo. 

Since we use the redshift of infall as interaction red- 
shift, mergers and flybys in our simulation will occur well 
before any observable encounter-induced features within 
the galaxies. For instance, the peak in the merger rate 
occurs for z ^ 5 (see F ig. SI) in stead of ~ 2 a s found 
elsewhere (|Madau et al . 1996; Rvan et all 120081 ). Note, 
though, that this definition of recording a merger as the 
redshift of infall is in line with the techniques used in the 
literature to construct mcrgcrtrees using FOF halos. 

3.1. Examples of halo evolution within the network 

In Fig. [1] we show the halo interaction network for two 
Milky- Way type halos of 10^^ Mg/h at redshift zercEI. 
These are two drastically different halos, even though 
they have the same mass a,t z — 0. The first halo formed 
at z — 9.3 and had a very active merger history; by 2; = 0, 
the assembly of this halo required 434 mergers. The sec- 
ond halo formed earlier, at z = 12.5, yet only required 
111 mergers - a 'quiet' evolution. Much of the mass in 
the quiet halo is assembled early through smooth accre- 
tion of dark matter; at z '-^ 5, the hectic halo has a mass 
of 2 X 10^° Mg/h while the quiet halo is 10 times more 
massive. The hectic primary halo itself directly survives 
77 mergers, while the quiet halo undergoes only 27. Such 
divergent pasts can leave an imprint on the structure and 
properties of the galaxies contained in these halos - in 
the star formation history, the morphology, and perhaps 
the central black hole mass. We will explore this topic in 
a separate paper (Sinha & HoUey-Bockelmann, in prep). 

® We used the total mass of the main halo, including all the 
subhalos, as proxy for the virial mass, and assumed a spherical 



3.2. Finding Flybys 

One major goal of this paper is to characterize flyby 
interactions - interactions that do not end with one halo 
accreting another. In principle, there are three classes of 
close flyby interactions: 

• Grazing ~ Two primary halos approach on con- 
vergent trajectories, interpenetrate for at least half 
a crossing time, and then separate as two distinct 
primary halos once more. 

• External - Same as above but the primary halos 
remain distinct at all times. This can also apply to 
halos that are at the same hierarchy levels within 
the same container halo (e.g., two subhalos of the 
same main halo). Since technically this can include 
every other halo in the volume, it is useful to define 
a maximum pericenter distance between two halos 
when classifying this type. 

• Internal - A halo at a higher hierarchy level passes 
close to the center of its containing halo, e.g., a 
sub-sub halo goes through the central regions of its 
containing subhalo. This is synonymous with the 
decay of satellite orbits after a merger. 
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Kazantzidis et al.i i2009f ) , so we will not discuss them 
here. External flybys are naturally distant encounters, 
where the separation between halo centers, Rsop, is 
larger than the virial radius of each halo. Since the 
perturbati on in the potential induced by a halo flyby 
is oc Rs(,p (|Vesperini fc Weinbergl l2000l ) . external flybys 
excite comparatively weak perturbations in the galactic 
potential. Since we are interested in potentially transfor- 
mative interactions, we will focus on grazing flybys. For 
the remainder of this paper, a 'flyby' means a grazing 
encounter. 

We require that the subhalo in flyby remain in the pri- 
mary halo for at least half a crossing time at the time 

of infall. The crossing time, icross — y^Rvir'^/GMvir, 
is independent of mass. Because we define halos by 
an overdensity, Avir = 200 x pcrit, Across simplifies to 
0.1 X H(z)~^. Our flyby definition is physically moti- 
vated in that we concentrate on longer-duration flybys 
that are more likely to leave an lasting imprint on the 
structure of the halos. However, this definition does ex- 
clude some rapid, transient events - hence, the results 
presented in this paper are conservative estimates of the 
true flyby rate in the Universe. 

For grazing encounters, the interaction will always be 
a flyby for very large relative velocities (urci >> Wcsc of 
the combined halo masses). However, we further com- 
pute the velocity dot product for the main halo centers 
in the center of mass frame at the snapshot where the ha- 
los are separate, (wi — Vfj_) ■ (v2 — w^)- If the dot product 
is negative, then the halos were approaching before they 
interpenetrate. We found that the flybys with negative 
V ■ V are slower and sink deeper into the main halo than 

halo to get the virial radius. 
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Figure 1. Left: Merger tree for a hectic halo. The mass of the halo at 2 = and z = 4.5 is 1.8 x 10^^ Mo/h and 2 x 10^° Mo/h 
respectively. Red represents a primary halo, blue represents a subhalo, green represents a subhalo that is going to dissolve inside 
the main halo, and black represents a halo that skips at least one snapshot. The size of the circle represents the halo's log-scaled 
mass with respect to the mass of the final main halo at z — 0. We only show those halos that are at least 1/10 of the primary 
mass. The corresponding redshifts are shown on the left. The fiybys are shown on the right with yellow circles; the size of the 
circle is scaled to show the mass ratio of the flyby. The main halo undergoes 77 mergers and 21 fiybys; however, we can only 
see a specific type of flyby in this mergertree - those relatively rare fiybys that eventually result in a merger. Right: Merger 
tree for a quiet MW mass halo. Fiybys are shown as yellow circles on the right. If a subhalo actually survives until z — 0, then 
it will not appear in this mergertree. At z — 0, the halo mass is is 10^^ Mq/Il , while at z = 4.5, it is 2 x 10^^ M0/h; this halo 
assembled most of its mass early via smooth accretion. 
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Figure 3. The solid white circle marks the virial radius of 
a typical halo extracted from the 512^ simulation (left) and 
the 1024^ simulation (right) at z; ~ 1.8. Dashed white cir- 
cles locate the other FOFs in the volume. Identical cubes of 
side 1 Mpc/h were extracted from both the simulations and 
projected to make this plot. This halo experienced the same 
number of merger and flybys throughout the run. 
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Figure 2. This figure shows a 5:1 flyby between a Milky Way 
type halo and a satellite halo. The flyby starts at z = 1.1 
and lasts until z — 0.3, for a duration of 4.7 Gyr. The gray 
and red points show the particle distribution of the FOF and 
the subhalo. The black and the dark red circles show the 
virial radius of the main halo and the subhalo throughout 
the encounter while the blue line shows the trajectory of the 
subhalo center of mass. This flgure also shows the artiflcial 
mass-loss experienced by subhalos as they pass close to the 
main halo center. 

the flybys with positive v • v. Fig. [2] shows an example 
of a deep flyby with a mass ratio of 5:1 from the testbed 
simulation. The event begins at z = 1.1 and completes 
by z = 0.3 - a duration of ~ 4.7 Gyr. At infall, the main 
halo and subhalo virial radii are 260 and 150 kpc/h 
respectively. We find the two halo centers are closest at 
z ^ 0.8 with a separation of ~ 100 kpc/h; the subhalo 
itself overlaps the main halo center here, which suggests 
that this encounter can perturb even the innermost re- 
gions of the main halo. Our technique is designed to 
identify this type of strongly-interacting, long-duration 
flyby, as it has the most potential to transform the over- 
all structure of each halo. 

4. RESULTS 

4.1. Mass and time resolution effects 

In this section, we estimate the effect of particle reso- 
lution and snapshot timing resolution on the flyby rate. 
We have used a 512"^ simulation and a 1024^^ simulation 
that is exactly the same phase for the results in this sec- 
tion. 

In order to determine the sensitivity of flyby rates 
to numerical resolution, we conducted two experiments 
that degraded our fiducial simulation. In the first ex- 
periment, we explored the effect of particle number on 
the frequency of flybys with a 512'^ simulation that was 
seeded with precisely the same initial random phases as 
the 1024'^ run (Roman Scoccimarro, private communi- 
cation) . Figure [3] centers on a typical MW halo in the 
volume - this halo experienced no difference in the num- 
ber of mergers or flybys over its entire history. In the 



entire volume, the difference between the number of fly- 
bys in the two simulations is less than 1 percent, when 
we consider the halos that are well-resolved in both simu- 
lations (i.e., down to the minimum halo mass in the 512'^ 
volume of M < 6.7 x lOi^Mg/h). In addition, the num- 
ber of flybys in two-dimensional mass and redshift bins 
are identical within Poisson error for 95% of the bins for 
these two simulations. So we find that particle resolution 
does not affect flyby rates. 

The second series of tests concerned the timing res- 
olution of the simulation. In the 512'^ simulation, we 
wrote snapshots in equal increments of log a, but never 
allowed more than 200 Myr between snapshots - this 
upper limit ensures that encounters for MW mass ob- 
jects were resolved over the course of the run. With 161 
snapshots, we made two new halo interaction network 
by skipping one and two snapshots respectively. Thus, 
we have three samples built from 161, 80 and 54 snap- 
shots respectively to quantify the effects of snapshot tim- 
ing resolution on the flyby rate. We analyze these three 
samples and compare the total number of flybys between 
halos. We find that to get convergence in the flyby rate, 
the snapshot outputs need to resolve the 'typical' flyby 
duration. Since a flyby is tagged as a main halo sub- 
halo main halo transition, and we require flybys to 
last at least half a crossing time, the snapshots must be 
able to resolve events on timescales of ^crossing- For the 
simulations with 161 and 80 snapshots, the crossing time 
was always resolved by snapshots and resulting number 
of flybys differ by less than 1% and are within Poisson 
errors. However, the 54 snapshot simulation did not re- 
solve the crossing time for most redshifts, and as a result 
the number of flybys is ^ 20% lower. 

As another test, we re-simulated the volume between 
z — 3 — 2 with a snapshot timing resolution of 25 Myrs, 
which yielded 45 snapshots. The crossing times at z = 3 
is ^ 330 Myrs, so it is very well resolved in this simula- 
tions. We take a similar approach as above, creating a 
series of progressively coarser halo interaction networks 
by skipping over snapshots. As before, we find that the 
number of flybys is consistent within Poisson error when 
we resolve the crossing time and drops dramatically for 
coarser simulations. Thus, our definition of flybys pro- 
duces consistent flyby rates as long as the snapshot out- 
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puts can resolve the crossing time. Otherwise, the flyby 
rates reduce as the time resolution degrades. This implies 
simulations with insufficient outputs will underestimate 
the true flyby rate. Assuming that the time between 
snapshots is at most a quarter of the crossing time (such 
that flybys last for at least 2 snapshots) , we find that the 
total number of snapshots required between redshifts 20 
and to be ~ 132. This is in agreement with recent 
findings that show at least 128 snapshots a re required to 
obtain a robust estimate of halo masses ([Benson et al.l 

[Mil. 

4.2. Frequency of flybys 

Now that we have a way to track and classify interac- 
tions between halos, it is useful to determine how com- 
mon close flyby encounters are compared to mergers. 
This is an important first step in determining whether 
current semi-analytic and high resolution N-body mod- 
els could be missing a significant number of galaxy- 
transforming interactions. Note that in this census, we 
are simply keeping a tally of dark matter halo mergers 
and close halo flybys, regardless of how much they may 
perturb the potential of any galaxy embedded within. In 
the next paper, we will concentrate on quantifying the 
perturbation caused by these types of encounters. We 
checked to see if there are multiple flybys between the 
same halo pairs and found that ~ 70% of the flybys do 
not recur. About 20% of flybys eventually become a 
merger, ~ 6% are repeat flybys. Thus, flybys primarily 
represent one-off events between two halos. 

In Fig. UJ we present the number of flybys and mergers 
per Gyr per cubic Mpc as a function of redshift. Mergers 
dominate the overall rate of interactions for 12 < z < 4 
by about an order of magnitude. However, at z '~ 14, 
we see tentative evidence of an increase in the relative 
number of flybys compared to mergers. At low redshift, 
2^3, flybys start to become more prevalent. The im- 
pact of this on observed merger rate estimates may be 
profound: in large-scale surveys, the rate of close-pair 
galaxies is often taken as a proxy for merger rate, but if 
we assume that 50% of these halo flybys appear as galaxy 
flybys, then Fig. 2] implies that at z < 3, such galaxy- 
pair counts are going to be contaminated by flyby-halos 
at least at a 20-30% level. Flybys also result in bimodal 
kinematic galaxy distribution near the virial radius, in 
that some galaxies are infalling for the flrst time, while 
other galaxies that were once much deeper inside the 
cluster potential are now leaving. Conceptually, these 
flyby galaxies are similar to 'backsplash galaxies' (e.g . , 
IGill et al.l[2005t IMoore et al.l[200l : IKnebe et al.ll2011bri . 

While Fig. |4] shows the global fiyby rate, we need to 
know the rate of interactions as a function of both halo 
mass and redshift to assess the impact flybys may have on 
halo evolution. Interaction rates can be expressed either 
as 'per object' or 'per volume'. Here, we plot the rates 
per halo to reduce the additional dependence of the halo 
mass function evolution. In Fig. [5] we show the number 
of flybys and mergers on a per halo per Gyr basis, while 
the ratio of flybys to mergers is seen in Fig. [51 We see 
that flybys are more frequent for higher-mass halos at all 
redshifts, consistent with the ACDM framework. From 
Fig. ini we can directly see that the number of flybys be- 
comes comparable or larger than the number of mergers 
for halos > 10^^ Mq /h for z < 2. Such halos are expected 
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Figure 4. The top panel shows the number of mergers (blue, 
dashed line) and flybys (red, solid line) per co-moving Mpc^ 
per Gyr as a function of redshift. The shaded region shows 
the Poisson error on the number of interactions. The bottom 
panel shows the ratio between the number of flybys to mergers 
vs redshift. Mergers dominate flybys in the redshift range ~ 
4—12 but flybys become relatively more common at z < 3. We 
can also see tentative evidence for an increase in the relative 
fraction of flybys compared to mergers at z > 14. However, 
this increase is also due to the low number of interactions and 
is borne out by the large Poisson error at those high redshifts. 
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Figure 5. This figure shows the number of mergers (left) and 
flybys (right) per halo per Gyr as a function of primary halo 
mass and redshift. The dashed line shows the mass accretion 
history of a typical Milky- W ay type obtained from our simu- 
lations (jWechsler et al.ir2002f ). The number of flybys increases 
with the primary halo mass for all redshifts, consistent with 
ACDM. For z < 3, halos above lO^^Mo/h are have flyby 
rates greater than 100 per Gyr. Such a high flyby rate (in 
addition to mergers) means these halos are unlikely to be in 
equilibrium. 

to host galaxies - and the effect of flybys should leave an 
imprint on the observable properties. We also find ten- 
tative evidence that flybys are comparable in number to 
mergers at the highest redshifts, z > 14 - however, the 
numbers are affected by Poisson error and we can not 
conclusively say that flybys dominate mergers at those 
highest redshifts. We are simulating a bigger boxsize to 
look at the relative importance of flybys for z > 14. 

In the preceding mergers versus flybys comparison fig- 
ures, we used all interactions irrespective of mass-ratio, 
q. Now, it is interesting to see how the frequency of in- 
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Figure 6. Ratio of number of flybys to mergers as a function 
of primary halo mass and redshift. As we saw in Fig. 31 
mergers dominate flybys by an order of magnitude for 12 < 
2 < 4. At lower redshifts, however, flybys start becoming 
more prevalent and by z ~ 2, flybys are at comparable or 
even larger for all halos above 10^^ Mq/Yi. 
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Figure 7. The average number of interactions per Gyr for a 
MW type halo (selected at z = 0) as a function of mass-ratio 
(q) and redshift. Mergers are plotted on the left panel while 
flybys are on the right panel. As seen in the previous plots, 
mergers dominate flybys for all q values in the redshift range 
12 — 2 but at flybys become comparable in number at lower 
redshifts. 

teractions compare as a function of both z and q for a 
typical MW halo. We select MW halos at z = using 
three criteria: a) halo mass between 1 and 2 x 10^^ M© /h 
and b) the halo has always been a primary halo and c) 
the main progenitor of the halo can be traced back to 
z > 6. We found 178 such halos in our simulation box. 
In Fig. [71 we show the average frequency of interactions 
per Gyr for these halos as a function of q and z. Mergers 
are dominant for 12 < z < 2 for all q; however, flybys 
take over for lower z. For these MW type halos, flybys 
and mergers happen with a similar frequency for z < 2 - 
a trend that we saw in Fig. |6l 

4.3. Fits to interaction rates 

Merger rates are usually calculated in two different 
ways: as a rate i) per parent halo ii) per child halo. These 
two rates physically mean slightly different things - the 
first one is tied to the probability of a halo of a given mass 
at some z having a merger. Observationally this is mea- 



sured by the merger fraction through pair counts. The 
second one is the probability that a halo had a merger 
in the past - observationally identified by morpholog- 
ical features. Since we track subhalos, we can distin- 
guish between the beginning of a merger (i.e., a main 
halo— i>subhalo transition - tied to the merger rate per 
parent halo) and the end of a merger (subhalo— ?>main 
halo - tied to the merger rate per child halo). These two 
methods of measuring the merger rate can yie ld some- 
what different results (see iHopkins et al.l 12010( 1 . In this 
section we will use the time of primary— >subhalo transi- 
tion as a merger time for the primary halo. 

To compare to merg er rates found i n the literature, we 
use the analytic fit of iWetzel et al.l ()2009l ) and classify 
interactions per halo per Gyr in the form: 



^events 
"■halodt 



= A(l + z)", 



(1) 



where, ricvonts is the number of mergers or flybys for a 
parent halo in some mass bin, nhaio is the number of halos 
in the same mass range and dt is the time interval (in 
Gyrs) between consecutive snapshots. We consider three 
different mass ranges for the parent halos, 10^" — 10^^, 
10" - 10^2 and 10^2 - IQia Mg/h. We limit the redshift 
range from z = 6 to 1.0. The fits for mergers and flybys 
are presented in Table [H Overall, the values of A for 
flybys and mergers are comparable, while a is higher for 
mergers. Correspondingly, the flyby rate increases more 
slowly than the merger rate with increasing z. 



Table 1 

Fits to Eqn.[T]for three main halo mass ranges. The redshift is 
noted when the subhalo first appears inside the main halo. The 
reduced x-squared values for the fits are also shown in the table. 



Mo 




10^1 - 10^2 


10^2 10^^ 




A 


1.33 X 10-3 


1.61 X 10-2 


3.84 X 10-1 


Mergers 


a 


2.40 


2.58 


1.16 






10.6 


3.05 


1.45 




A 


1.08 X 10^3 


1.92 X 10-2 


2.59 X IQ-i 


Flybys 


a 


1.02 


0.75 


0.68 






2.33 


1.63 


0.88 



5. DISCUSSION 

We present a technique for tracking the dynamical his- 
tory of halos in an N-body simulation. Since many of 
these interactions are close flybys and not mergers, we 
call this ensemble the halo interaction network. We find 
that most of the anomalies common to merger tree con- 
struction arise from limitations in the halo finding algo- 
rithm, and we correct for those in creating the halo in- 
teraction network. For example, density contrast based 
methods like SUBFIND are unable to track a subhalo 
when it reaches a denser central region. In most cases, 
such close pericenter passages manifest as severe mass 
loss in the subhalo (with a fraction of the mass regained 
when the subhalo is farther away from the host center), 
and at times, the entire subhalo falls below the mini- 
mum particle limit and is not identified at all. There are 
also issues of nomenclature - in a near equal-mass major 
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merger, it is unclear which halo is the dominant one, i.e., 
which halo is to be labeled the 'main' halo. The algo- 
rithm we have presented here tracks several pathological 
cases and corrects them. We find that even after fixing 
all such cases, there still remain halos that are parentless 
or childless at a few-10% level for halos comprised of < 
80 particles. Even with very finely spaced snapshot out- 
puts, we can not eliminate all of these ill-behaved halos. 
Based on our findings, we recommend that semi-analytic 
studies based on N-body mergertrees use only those ha- 
los with at least 3 — 4 times the formal particle limit. In 
our simulations, we formally require 20 particles to iden- 
tify a halo, but all the results that have been reported 
here involve halos with at least 100 particles. 

For robust halos, our mergertree technique resolves 
anomalies related to halo definition and problems with 
halo identification; however, the code can still be im- 
proved using maximal bi partite graph matching (e.g., 
iHopcroft fc Richardlll973[ ). Given two sets of halos at 
two different snapshots and an associated cost (say, the 
inverse of the binding energy rank) for matching any two 
of them, the maximal matching algorithms can be used 
to determine a match. Such a technique would produce 
the least number of anomalies and spurious initial as- 
signments in a mergertree, and we will implement it in a 
future version of the mergertree code. 

One major thrust of this paper is to include fiyby in- 
teractions in our halo interaction network. In this paper, 
we present the first census of halo flybys across a wide 
range of halo mass and redshift. We classify a fiyby us- 
ing spatial and dynamical information. Conceptually, a 
grazing fiyby occurs when a halo undergoes the transi- 
tion from a main halo — >-subhalo— >■ a main halo and is 
fundamentally different from mergers. In a simulation, it 
is relatively straightforward to identify such chains since 
the current kinematics and future behavior of a given 
interaction is fully determined. We note that because 
flybys imply a population of galaxies that were once well 
within the virial radius of the main halo but are now 
outside it, they may be related to so called 'backsplash' 
galaxies. The morphological and kinematical features of 
these galaxies are likely to be different compar ed to other 
galax i es that are infalling for the first time ()Gill et al.l 
I2OOI IMoore"eFan[20Ql IKnebe et al.ll2011bD . In some 
cases these flybys can even result in subhalos switch- 
ing their main halo an d appear as 'renegade' subha- 
los (IKnebe et al.l l20lTal ). While there has been a sig- 
nificant amount of work related to ' backsplash' galaxies 
in clusters (e.g. . 
Gill et all 120051:' 



Balogh et al 



Moore et alj 



[2OOOI: iSanchis et al 



2004t iMamon et al 



2002 



20041: 



Knebe et al.ll2011bl: iMahaian et al.il201lD . unfortunatelv 
there are no systematic studies for lower mass halos. 

We suspect fiyby interactions have been unappreci- 
ated in both observational and theoretical studies that 
attribute morphological/structural changes to mergers. 
For example, the pre-dominance of low-mass bulge- 
dominated galaxies (Mstdiar ^ 10^°MQ/h) may be po- 
tentially explained by perturbations caused by flybys, 
since low-mass gala xies do not suffer eno ugh mergers to 
explain the bulges (jWeinzirl et al.l l2009l ). Indeed, fly- 
bys can cause a fractional change in the binding energy 
of the pr imary halo, AE/E > 1% - this may excite 
bars (see iBer entzen et al.l 12004 and references therein) 
and drive gas into the central regions of the galaxy, cre- 



ating a pseudo-bulge. Observationally, barred and un- 
barred galaxies can not be separated on some inherent 
galaxy property (other than the actual presence/absence 
of a bar); this implies that the formation or destruc- 
tion of a bar is li kely lin ked with some external process, 
such as a flyby (jSellwo od 2010). In addition to form- 
ing bars/pseudo-bulges, flybys could create density en- 
hancements via shocks in the hot halo gas of the primary 
galaxy. Such density enhancements will increase the 
cooling rates in the halo by ^ an order of magnitude (as- 
suming bremsstrahlung and a factor of 4 increase in den- 
sity). In principle t his can cause the hot gas to cool an d 
replenish the disk (jSinha fc HoUev-Bockelmannl I2009D . 
In a similar spirit, flyby-induced star formation has been 
recently invoked to explain the steepness of the [a/Fe]-a 
relation (jCalura fc Mcnci 20ill). Flybys can also cause 
various other morphological features or transformat ions , 
e.g., spiral to SO galaxy in gr oups (Bekki fc Cou ch,.20Tll) . 
spin flips in the inner halo (IBett fc Frenkll2012D or spi- 
ral a rms in the galactic disk fe.g.. ITutukov fc Fedoroval 
I2OO6I) . 

6. CONCLUSIONS 

In this paper we report on a new class of interactions 
~ flybys, that occur frequently. Most of these flybys are 
one-off events - one halo delves within the virial radius 
of another main halo, separates at a later time and does 
not return. From our testbed simulation, we can see 
a hint that flybys dominate over mergers at very high 
redshift (z ^ 14) and have relatively shallower depen- 
dence on z compared to mergers (see Table [T|). We find 
that most fiybys are one-off events and about 70% of 
the fiybys do not ever return. Flybys are then, a largely- 
ignored type of interaction that can potentially transform 
galaxies. Unfortunately, most semi-analytic methods of 
galaxy formation are designed only to use mergers and 
thus can not account for the effects of flybys directly. In 
principle, a flyby can create long-lasting features to ap- 
pear in the secondary halo and affect the evolution of 
that halo. In general, slow flybys cause a larger pertur- 
bation compared to a fast one, and such features can 
persis t even when the pertu rbing halo has moved far 
away (iVesperini fc Weinberd '2000). If flybys cause per- 
ceptible changes in galactic structure, then those devi- 
ations should manifest in the oldest galaxies. Upcom- 
ing space telescopes like JWST will target high redshift 
(z ^ 7) galaxies to unravel clues for galaxy formation 
and may be able to detect flyby signatures. 

In conducting our flyby census, our testbed simulation 
did not have adequate resolution to capture the per- 
turbative effect of a flyby on the dynamics of the halo 
central region. We are currently exploring the effect 
of flyby interactions on galaxy structure and morphol- 
ogy using both high resolution N-body simulations and 
semi-analytic estimates (Sinha fc HoUey-Bockelmann, in 
prep). 

This work was conducted in part using the resources 
of the Advanced Computing Center for Research and 
Education at Vanderbilt University, Nashville, TN. We 
also acknowledge support from the NSF Career award 
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Pleiades cluster at NASA Advanced Supercomputing Di- 
vision. We would like to thank the anonymous referee for 
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also like to thank H.R. Madhusudan for pointing out bi- 
partite graph matching algorithms. 



REFERENCES 



Balogh, M. L., Navarro, J. F., & Morris, S. L. 2000, ApJ, 540, 113 

Barnes, J. E. 2002, MNRAS, 333, 481 

Barnes, J. E., & Hernquist, L. 1992, ARA&A, 30, 705 

— . 1996, ApJ, 471, 115 

Behroozi, P. S., Conroy, C, & Wechsler, R. H. 2010, ApJ, 717, 

379 

Behroozi, P. S., Wechsler, R. H., & Wu, H.-Y. 2011, ArXiv 

c-prints 

Bckki, K. 2008, MNRAS, 388, LIO 

Bekki, K., & Couch, W. J. 2011, MNRAS, 415, 1783 

Benson, A. J., Borgani, S., De Lucia, G., Boylan-Kolchin, M., & 

Monaco, P. 2011, MNRAS, 1924 
Berentzen, I., Athanassoula, E., Heller, C. H., & Fricke, K. J. 

2004, MNRAS, 347, 220 
Bett, P. E., & Frciik, C. S. 2012, MNRAS, 420, 3324 
Boylan-Kolchin, M., Springel, V., White, S. D. M., Jenkins, A., & 

Lemson, G. 2009, MNRAS, 398, 1150 
Bridge, C. R., Appleton, P. N., ConseUce, C. J., Choi, P. I., 

Armus, L., Fadda, D., Lainc, S., Marleau, F. R., Carlberg, 

R. C, Hclou, C, & Yan, L. 2007, ApJ, 659, 931 
Calura, F., & Mcnci, N. 2011, MNRAS, 413, LI 
Casasola, V., Bettoni, D., & Galletta, G. 2004, A&A, 422, 941 
Covington, M., Dekel, A.. Cox. T. J., Jonsson, P., & Primack, 

J. R. 2008, MNRAS, 384, 94 
Cox, T. J., Di Mattco, T., Hernquist, L., Hopkins, P. F., 

Robertson, B., & Springel, V. 2006, ApJ, 643, 692 
Cox, T. J., Jonsson, P., Somerville, R. S., Primack, J. R., & 

Dekel, A. 2008, MNRAS, 384, 386 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, 

ApJ, 292, 371 
De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 
Dubinski, J., Mihos, J. C, & Hernquist, L. 1999, ApJ, 526, 607 
Genel, S., Genzel, R., Bouche, N., Naab, T., & Sternberg, A. 

2009, ApJ, 701, 2002 
Genel, S., Genzel, R., Bouche, N., Sternberg, A., Naab, T., 

Schreiber, N. M. F., Shapiro, K. L., Tacconi, L. J., Lutz, D., 

Cresci, G., Buschkamp, P., Davies, R. I., & Hicks, E. K. S. 

2008, ApJ, 688, 789 
Gill, S. P. D., Knebe, A., & Gibson, B. K. 2005, MNRAS, 356, 

1327 

Guo, Q., & White, S. D. M. 2008, MNRAS, 384, 2 
Hernquist, L., & Mihos, J. C. 1995, ApJ, 448, 41 
HoUey-Bockelmann, K., & Richstone, D. O. 2000, ApJ, 531, 232 
Holmberg, E. 1941, ApJ, 94, 385 

Hopcroft, J. E., & Richard, M. K. 1973, SIAM Journal on 

Computing, 2, 225 
Hopkins, P. F., Croton, D., Bundy, K., Khochfar, S., van den 

Bosch, F., Somerville, R. S., Wetzel, A., Keres, D., Hernquist, 

L., Stewart, K., Younger, J. D., Genel, S., & Ma, C. 2010, ApJ, 

724, 915 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Martini, 
P., Robertson, B., & Springel, V. 2005, ApJ, 630, 705 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 
Robertson, B., & Springel, V. 2006, ApJS, 163, 1 

Johnston, K. V. 1998, ApJ, 495, 297 

Kazantzidis, S., Zentner, A. R., Kravtsov, A. V., Bullock, J. S., & 

Debattista, V. P. 2009, ApJ, 700, 1896 
Klypin, A., Gottlober, S., Kravtsov, A. V., & Khokhlov, A. M. 

1999, ApJ, 516, 530 
Knebe, A., Libeskind, N. L, Doumler, T., Yepes, G., Gottlober, 

S., & Hoflfman, Y. 2011a, MNRAS, 417, L56 
Knebe, A., Libeskind, N. L, Knollmann, S. R., Martinez- Vaquero, 

L. A., Yepes, G., Gottlober, S., & Hoffman, Y. 2011b, MNRAS, 

412, 529 

Knollmann, S. R., & Knebe, A. 2009, ApJS, 182, 608 



Komatsu, E., Dunklcy, J., Nolta, M. R., Bennett, C. L., Gold, B., 
Hinshaw, G., Jarosik, N., Larson, D., Limon, M., Page, L., 
Spergel, D. N., Halpern, M., Hill, R. S., Kogut, A., Meyer, 
S. S., Tucker, G. S., Wciland, J. L., WoUack, E., & Wright, 

E. L. 2009, ApJS, 180, 330 

Lacey, C, & Cole, S. 1993, MNRAS, 262, 627 

Laine, S., van der Marel, R. P., Rossa, J., Hibbard, J. E., Mihos, 

J. C, Boker, T., & Zabludoff, A. L 2003, AJ, 126, 2717 
Lin, L., Patton, D. R., Koo, D. C, Casteels, K., Conselice, C. J., 

Faber, S. M., Lotz, J., Willmer, C. N. A., Hsieh, B. C, Chiueh, 

T., Newman, J. A., NovaJc, G. S., Weiner, B. J., & Cooper, 

M. C. 2008, ApJ, 681, 232 
Lotz, J. M., Jonsson, P., Cox, T. J., & Primack, J. R. 2010, 

MNRAS, 404, 575 
Ludlow, A. D., Navarro, J. F., Springel, V., Jenkins, A., Frenk, 

C. S., & Helmi, A. 2009, ApJ, 692, 931 
Maciejewski, M., Colombi, S., Springel, V., Alard, C, & Bouchet, 

F. R. 2009, MNRAS, 396, 1329 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., 
Steidel, C. C, & Fruchter, A. 1996, MNRAS, 283, 1388 

Mahajan, S., Mamon, G. A., & Raychaudhury, S. 2011, MNRAS, 
416, 2882 

Mailer, A. H., Katz, N., Keres, D., Dave, R., & Weinberg, D. H. 

2006, ApJ, 647, 763 

Mamon, G. A., Sanchis, T., Salvador-Sole, E., & Solanes, J. M. 

2004, A&A, 414, 445 
Micic, M.. HoUey-Bockelmann, K., & Sigurdsson, S. 2011, 

MNRAS, 450 

Micic, M., Holley-Bockelmann, K., Sigurdsson, S., & Abel, T. 

2007, MNRAS, 380, 1533 

Mihos, J. C, & Hernquist, L. 1994, ApJ, 431, L9 

— . 1996, ApJ, 464, 641 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Moore, B., Diemand, J., & Stadel, J. 2004, in lAU CoUoq. 195; 

Outskirts of Galaxy Clusters: Intense Life in the Suburbs, ed. 

A. Diaferio, 513-518 
Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A. 1996, 

Nature, 379, 613 
Muldrew, S. I., Pearce, F. R., & Power, C. 2010, MNRAS, 1540 
Naab, T., Jesseit, R., & Burkert, A. 2006a, MNRAS, 372, 839 
Naab, T., Khochfar, S., & Burkert, A. 2006b, ApJ, 636, L81 
Neyrinck, M. C, Gnedin, N. Y., & Hamilton, A. J. S. 2005, 

MNRAS, 356, 1222 
Okamoto, T., & Habe, A. 2000, PASJ, 52, 457 
Purcell, C. W., Kazantzidis, S., & Bullock, J. S. 2009, ApJ, 694, 

L98 

Quinn, P. J., Hernquist, L., & FuUagar, D. P. 1993, ApJ, 403, 74 
Roberts, T. P., Warwick, R. S., Ward, M. J., & Murray, S. S. 

2002, MNRAS, 337, 677 
Robertson, B., Bullock, J. S., Cox, T. J., Di Mattco, T., 

Hernquist, L., Springel, V., & Yoshida, N. 2006a, ApJ, 645, 986 
Robertson, B., Cox, T. J., Hernquist, L., Franx, M.. Hopkins, 

P. F., Martini, P., & Springel, V. 2006b, ApJ, 641, 21 
Ryan, Jr., R. E., Cohen, S. H., Windhorst, R. A., & Silk, J. 2008, 

ApJ, 678, 751 

Sales, L. V., Navarro, J. F., Abadi, M. G., & Steinmetz, M. 2007, 

MNRAS, 379, 1475 
Sanchis, T., Solanes, J. M., Salvador-Sole, E., Fouque, P., & 

Manrique, A. 2002, ApJ, 580, 164 
Schweizer, F. 1986, Science, 231, 227 
Sellwood, J. A. 2010, ArXiv e-prints 

Sellwood, J. A., Nelson, R. W., & Tremaine, S. 1998, ApJ, 506, 

590 

Sharma, S., & Johnston, K. V. 2009, ApJ, 703, 1061 

Simha, V., Weinberg, D. H., Dave, R., Gnedin, O. Y., Katz, N., 

& Keres, D. 2009, MNRAS, 399, 650 
Sinha, M., & Holley-Bockelmann, K. 2009, MNRAS, 397, 190 
Springel, V. 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005a, ApJ, 620, L79 
— . 2005b, MNRAS, 361, 776 

Springel, V., & White, S. D. M. 1999, MNRAS, 307, 162 
Springel, V., White, S. D. M., Jenkins, A., Frenk, C. S., Yoshida, 
N., Gao, L., Navarro, J., Thacker, R., Croton, D., Helly, J., 
Peacock, J. A., Cole, S., Thomas, P., Couchman, H., Evrard, 
A., Colberg, J., & Pearce, F. 2005c, Nature, 435, 629 
Springel, V., White, S. D. M., Tormen, G., & Kauffmann, G. 
2001a, MNRAS, 328, 726 



10 



Springel, V., Yoshida, N., & White, S. D. M. 2001b, New 
Astronomy, 6, 79 

Tacconi, L. J., Genzel, R., Smail, I., Neri, R., Chapman, S. C, 
Ivison, R. J., Blain, A., Cox, P., Omont, A., Bertoldi, F., 
Greve, T., Forster Schreiber, N. M., Genel, S., Lutz, D., 
Swinbank, A. M., Shapley, A. E., Erb, D. K., Cimatti, A., 
Daddi, E., & Baker, A. J. 2008, ApJ, 680, 246 

Taylor, J. E., & Babul, A. 2001, ApJ, 559, 716 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 

Toth, G., & Ostriker, J. P. 1992, ApJ, 389, 5 

Tutukov, A. v., & Fedorova, A. V. 2006, Astronomy Reports, 50, 
785 

Tweed, D., Devriendt, J., Blaizot, J., Colombi, S., & Slyz, A. 
2009, A&A, 506, 647 



Urrutia, T., Lacy, M., & Becker, R. H. 2008, ApJ, 674, 80 

van Dokkum, P. G. 2005, AJ, 130, 2647 

Vesperini, E., & Weinberg, M. D. 2000, ApJ, 534, 598 

Walker, I. R., Mihos, J. C, & Hernquist, L. 1996, ApJ, 460, 121 

Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V., 

&; Dekel, A. 2002, ApJ, 568, 52 
Weinberg, M. D., & Blitz, L. 2006, ApJ, 641, L33 
Weinzirl, T., Jogee, S., Khochfar, S., Burkert, A., & Kormendy, 

J. 2009, ApJ, 696, 411 
Wetzel, A. R., Cohn, J. D., & White, M. 2009, MNRAS, 395, 1376 
Younger, J. D., Hopkins, P. F., Cox, T. J., & Hernquist, L. 2008, 

ApJ, 686, 815 



APPENDIX 

HALO GENEALOGY: FINDING PROPER PARENT-CHILDREN PAIRS 

Our first step is to track the main halos through time. Since even simple isolated halos change mass and position in 
each snapshot, this is not a trivial task. Overall, we link main halos and subhalos with different techniques. In most 
cases, we can assign the proper parent for a main halo in the following way: 

• Gather all the particles in all main halos at two consecutive snapshots. Note that this includes the particles in 
the subhalos because some of these particles may have joined (or left) the main halo between these two redshifts. 

• Working backwards from a later snapshot, we compare the particle ids in one halo with the ids in all primary 
halos at the previous snapshot. For computational efficiency, once a particle id is matched across two snapshots, 
we store the actual subhalos at the particle level. 

• Assign the parenlQhalo as the one with the largest number of common particles. 

• Repeat for all remaining primary halos at the later snapshot 

• If several child halos claim the same parent (i.e. the parent halo may have split apart), the parent is assigned 
to the most massive child halo. In other words, a parent may have at most one child halo in our scheme. In 
addition, one child halo can have at most one parent at this stage - though this may change at the subhalo 
matching step. 

Using this procedure, we assign ^ 90% of the main halos to a parent. If the halo is still unassigned, we treat it as 
though it were a subhalo. Subhalos are matched using a binding energy technique described below. 

In the SUBFIND algorithm, a subhalo is a self-bound group of particles. Therefore, an acceptable parent-child pair 
should have most of the hi ghly bound particles in common. We track the subhalo core (similar to Okamoto & Hab^ 
HOOO; Springel et al."2005c': ^ Lucia fc B laizot' 20071: iBovlan-Kolchin et al.l l2009') using a binding energy rank similar 
to .Bovlan-Kolchin et al.. (,2009ir For every particle that is common between two halos, we compute the quality - Q, 
as: 

2 = E ^^P-+ E ^.--xt , (Al) 

i<100 jXlOO 

where 'prev' and 'next' refer to the halo at the previous and later snapshot and TZ is the binding energy rank (1 for 
the most bound, 2 for the next and so on) of the common particle; this choice weights the core of a halo more strongly 
than the outer regions. We limit the number of particles in Q to 100, because in cases of severe tidal stripping, only 
the core remains - the rest of the subhalo is spread to the main halo background. If no child subhalo is found, we 
search two later snapshots using this binding energy rank method. This yields a parent-child match for over 99% of 
the halos. 

After the parent-child matching is done for all halos and subhalos, we manually check the assignments for validity 
and make adjustments if necessary. We found that the parent-child pair mis-assignments typically fell into two classes 
- those caused by major mergers, and those caused by numerical resolution. Major mergers caused < 1% of the 
mis-assignments while particle resolution resulted in the remainder of the poor parent-child pairs. During a major 
merger, the process may be so violent that a large parent subhalo is the only thing that survives relatively intact, 
while the main parent halo is mixed throughout the remnant. In this case, the parent halo and subhalo effectively 
'swap' categories after the merger. However, because our halo matching takes place before subhalo matching, our 
algorithm inaccurately assigns both the parent halo and the subhalo to the same main child halo after the merger. 
The consequence of this is that it eliminates an entire branch of the merger tree. To correct the pairing, we search 
for this pattern in the parent-halo pairs, demote the main parent halo to a child subhalo, and reinstate the missing 
branch. 



° Here a parent halo lies in the past. Note that Extended Press- 
Schechter usually has the opposite convention. 
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Figure 8. Schematic for finding the parents of subhalos passing through a dense halo center. Time increases vertically down- 
wards. As the subhalo passes through the central region, the density contrast becomes too low for SUBFIND to identify it; 
consequently, original particles from the subhalo blend with the halo background. At some later snapshot, when the subhalo is 
sufficiently far from the center, SUBFIND identifies it as a subhalo (left). Left alone, the parent-matching algorithm will lose 
the subhalo in step 1 to the main halo in step 2, and mistakenly identify a new parentless subhalo in step 3. We remedy this 
by checking that the subhalo in step 3 can be found in the same halo in an earlier snapshot. 

Tracking the network of interactions 

The parent-finding step described in the previous section simply links one halo to another between two timesteps, 
where possible. To make the halo interaction network, we must both uniquely track every halo from its formation to 
its destruction (or to the present day if it survives) and categorize every dynamical interaction between two or more 
halos. The procedure for tracking the halo is quite similar to our technique for pairing an individual parent and child 
halo, and amounts to constructing a merger tree; for completeness sake, however, we describe how we follow a halo 
through its lifespan. First, for every child halo at a given snapshot, we identify all the parents that comprise the child 
- these parent halos are all direct parents, but may each have joined the child halo at different redshifts. We choose 
the most massive parent as the primary halo, and assign unique halo IDs that are preserved as we step backwards in 
time along a particular halo branch. 

While the halo-pairing step addressed most anomalies, there is one anomaly that can only be addressed during 
this step: when a subhalo passes through a dense region of the halo, it seems to disappear because SUBFIND cannot 
differentiate it from the background halo. This means that when the subhalo reappears several snapshots later, it will 
seem brand new. To graft these apparently new subhalo branches to their proper place in the network, we first locate 
all the subhalos that do not have parents, and search among all subhalos up to five snapshots earlier for a suitable 
parent based on particle IDs. A potential match contains more than 50% of the particles of the orphaned subhalo and 
is not a primary parent of any other halo. Once the parent subhalo match is found, we attach the apparently orphaned 
subhalo branch to this parent subhalo (see Fig [8] for a schematic). 

These orphans occur in a very small mass range - so fixing th em does not dram atically increase the overall 99% 
successful of parent-child assignments. The anomalies described in ([Tweed et al.ll2009| ) are also present in our technique: 

• No child: None of the halo particles were found in any other halo, up to 3 snapshots into the future. 

• No parent: Either none of the halo particles were found in another halo up to 5 snapshots in the past, or the 
potential parent halo has already been matched with a child halo. 

• Transients: The halo appears to be a transient phenomenon - none of the halo particles are present in any future 
or past halo. We determined this to be largely an effect of numerical noise, as we will discuss in the next section. 

Halo tracking and particle number 

In Fig. iniwe show the cumulative distribution, over all snapshots, for the total number of parentless, childless and 
transient halos versus halo mass. To generate the plot, we bin in mass for all the snapshots in a specified redshift 
range and count the total number of halos that are either parentless or childless. This is divided by the total number 
of halos in the same range to give a fractional occurrence of the anomalies. 

It is clear that the number of childless and transient halos drop off very rapidly with particle number, and is essentially 
zero for halos resolved with > 80 particles, or about 4 times our nominal halo detection threshold. This implies that 
the typical particle number threshold 20 particles) for resolving halos is much too small to produce robust results. 
A different halo finder and mergertree technique has found similar results (Dylan Tweed, private communication). 
With these results in mind, we find that a halo is unlikely to be noise once it is above ~ 4 times the canonical particle 
number threshold for halo detection. For the science results in the paper we will employ an even more conservative 
limit of 100 particles per halo, 5 times our halo detection threshold. 
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Figure 9. The cumulative fraction of halo anomalies as a function of halo mass. The different lines represent parentless halo 
(thick solid), parentless subhalos (thin solid), childless (dashed) and transient (dash-dotted). The number of childless and 
transient halos drop off very sharply with mass, suggesting that both are caused by numerical noise. The parentless halos show 
a peak near our halo detection threshold - these are all the newly-formed halos and hence do not have a parent. The parentless 
subhalos show a very gradual decline with mass, and can be interpreted as a combination of two effects: a) timing resolution 
- a new halo forms and is accreted by another halo to become a subhalo between two snapshots b) a subhalo splits into two. 
The upper X-axis shows the number of resolution elements in a given halo mass; the cyan dotted line shows our recommended 
minimum halo particle threshold to limit transient behavior at or below the 1% level for reasonable timing resolution. The black 
dotted line shows the halo detection limit (20 particles), while the cyan line is the 4 times that value. 
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Figure 10. This plot shows the dependence of parentless halos for the z range 6 — 5, as a function of time resolution of 
snapshots. The red line shows our fiducial simulation while the blue shows the a re-simulation with 10 times as many snapshots. 
The number of parentless halos in the fine simulation reduces by factors of 4 — 5 compared to the coarse one. While the number 
of parentless subhalos is also reduced in the fine simulation, the drop is smaller. With such a high time resolution, we can reduce 
the fraction of parentless (sub)halos to the few % level for halos that are only twice the detection limit. To achieve the fine time 
resolution shown here, a simulation would have to output about 2000 snapshots for z = 20 to (assuming steps in log a). The 
transients/childless halos show virtually no difference between the two simulations and have been omitted from the plot. 



Halo tracking and time resolution 

Both parentless halos and parentless subhalos can occur; we expect most of the parentless halos to be legitimate, 
newly formed halos. However, some of the parentless subhalos may actually be caused by a newly formed halo that 
is quickly accreted by another primary halo in between two snapshot outputs. To test this hypothesis, we increased 
the output resolution of the simulation by an order of magnitude between z = 6 — 5, resulting in a snapshot every few 
timesteps. 

Fig. [TOl shows the percentage of halos that are parentless/childless in the redshift range z = 6 — 5 for two simulations 
with differing snapshot output resolution. Our fiducial simulation, which we label as 'coarse' in this plot, already has 
good timing resolution (see Section [5]) but still produces anomalies at the ^ 5 — 10% level in the low mass halos. 
The 'fine' simulation increases the snapshot output by a factor of 10 which reduces the level of anomalies by a factor 
of 2 — 3 over the entire mass range. Fig [10] shows that the overall number of parentless subhalos dropped with this 
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increased output frequency, though never to the level of parentless halos. The slope for parentless subhalos are similar 
in the two simulations; in fact we note that the overall shape of the parentless subhalos mimics the cumulative mass 
function, albeit with a much lower amplitude. We expect that increasing the time resolution even further will reduce 
the amplitude, but a fraction of parentless subhalos will still remain. In the 'fine' simulation, the fraction of parentless 
subhalos reduces to the 1% at only twice the halo detection limit, but such accuracy in the halo interaction network 
comes at the price of disk output - a full-scale cosmological simulation would need ^ 2000 snapshots (assuming log a 
output spacing) to achieve the timing resolution of the 'fine' simulation. The number of childless/transient halos is 
indistinguishable between the two simulations - showing that they are indeed dominated by numerical noise. 

Comparison of our technique to previous numerical work 

A mergertree simply identifies the links between components of a halo across redshifts. For particle-based codes, 
this means tracking halo particles by particle id across various snapshots. Nevertheless, the complexity of halo 
histories makes the task non-trivi al. As w e noted earlier, most halo-finders use some form of a density con- 
trast (e.g., Springcl ct al. 2001a; K nollmann & Kncbc 2009; Klvoin et al. 199ij) to group particles int o a subhalo 
(sec VOBOZ'(Nevrinck_^aL|(2005), EnLink (Sharma & Johnston 2009), HSF (Macieie wski et all 120091 ) and Rock- 
star ( Behrpozi et al- 20 111) for different approaches). Such a technique fails to identify subhalos close to the center of 
the host halo. Hence it is important to isolate such cases and re-assign the subhalos that 'disappear'. In this section, 
we compare the various strategies that have been employed to construct a merge rtree based on subhalos identified in 
N-bod y simulations. In particular, w e cor npare our mergertre e methods to those in Wechs ler et al.l (|2002D . IWetzel et al.l 
([20091, 'Bovlan-Kolchin ct alj ([20090 and lTweed et al.l (|2009l ). 

In Wcchslcr et al. (2002) a parent-child halo was assigned such that more than 50% of the particles in the parent end 
up in the child along with those 50% of particles constituting more than 70% of the particles in the parent that end 
up in any halo. This technique is also labeled as the 'most massive progenitor' method. Such an algorithm is prone to 
mis-assignment of parent-child halos in cases where a subhalo undergoes massive stripping or completely disappears 
inside a h ost halo . 

In Wetz el et al.l (|2009| ). the authors use the 20 most bound particles of a subhalo (that also occur in the potential 
child) to construct a parent-child pair. They compute the potential of these common particles based on all the particles 
in the parent and weight against large increases in mass or scale-factor. However, the weighting function is linear in 
mass (Eqn. 2) while the potential itself is squared (Eqn. 1). If a subhalo appears to have dissolved into the host halo, 
the large potential term could outweigh the limiting weighting function. As such, it is not obvious that this formulation 
would always pick the 'right' subhalo as a child. In addition, since N-body simulations have discrete outputs, weighting 
by the ratio of the scale-factors of the parent-child pair could change depending on the output frequency of the same 
si mulation. 

iTweed et al.l (|2009[ ) specifically made the mergertree to account for the cases we have discussed in this paper. They 
assign parent-child pairs by choosing the halo pairs with the maximum common mass as the main parent. In addition 
they identify (and remedy, if possible) three classes of anomalies: 1) subhalos without a parent 2) subhalo that becomes 
the parent of a primary halo and 3) subhalo and host halo swap identities. As we mentioned, we find two of these cases 
in our technique: anomaly 1 occurs under two circumstances - when the time resolution of the outputs is too coarse or 
when the subhalo appeared to have 'dissolved' in the host halo at some earlier snapshot. We correct for the dissolved 
halos in our interaction network. Anomaly 2 does not appear in our simulations since we first match main halos across 
two snapshots before matching subhalos. An analog of anomaly 3 occurs in our technique for major mergers - when 
SUBFIND switches the main halo and dominant subhalo definition across two snapshots. We have outlined our m ethod 
to fix these cases in Section lA.ll Overall, the approach described here and those used by iTweed et al.l (l2009[ ) agree 
qualitatively. We reiterate the cautionary statement from ITweed et all ([2009) that careful analysis has to be done to 
grow an 'accurate' mergertree. Although the anomalies that we have described and fixed represent < 5% of all halos, 
it is not clear to us what the cumulative effects of such anomalies are on semi- analytic approaches. 

Specifically the advantages in our mergertree technique are as follows: 

First: It uses different strategies to locate descendants for main halos and subhalos. Descendants of main halos are 
assigned based on the highest num ber of common particles, subhalo descendants are assigned according to the 
binding energy rank (see Eqn. [AT|1 . The rank computation is symmetric, i.e., a highly bound common particle 
always adds to the total rank. Using only common particles or binding energy rank leads to mis-assignments, 
e.g., a highly stripped subhalo will be lost even when the core survives or a main halo will get assigned to a 
subhalo that happens to contains some of the most highly bound particles from the previous step. 

Second: It carefully corrects the following pathological cases to the 0.1% level: 

-: a subhalo and main halo switching definition across timesteps 

-: a subhalo appearing to dissolve into the main halo due to massive tidal stripping (even though the core of the 
subhalo survives) 

-: when a subhalo (including the core) seems to disappear for multiple snapshots while passing close to the center of 
a main halo 



